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An overview is given of recent advances in the nonequilibrium statistical mechanics of quantum 
systems and, especially, of time-reversal symmetry relations that have been discovered in this con- 
text. The systems considered are driven out of equilibriuiu by time-dependent forces or by coupling 
•/j ' to large reservoirs of particles and energy. The symmetry relations are established for the exchange 

Cn , of energy and particles between the subsystem and its environment. These results have impor- 

tant consequences. In particular, generalizations of the Kubo formula and the Casimir-Onsager 
reciprocity relations can be deduced beyond linear response properties. Applications to electron 
' m4 , quantum transport in mesoscopic semiconducting circuits are discussed. 

^ ! I. INTRODUCTION 

Microreversibility is the symmetry under time reversal of the equations of motion for the microscopic particles 
j^ ' composing matter. This symmetry is a property of the electromagnetic interaction [l|. The motion is ruled by 
Newton's equations for the positions and velocities of the particles in classical mechanics and by Schrodinger's equation 
for the wavefunction in quantum mechanics. Classical mechanics emerges out of quantum mechanics in the limit where 
''rt de Broglie's wavelength becomes smaller than the spatial scale of features in the energy landscape where each particle 



mv 



where h is Planck's constant. Hence, the wavelength is shorter for heavier particles in systems at higher temperature. 
At room temperature, the wavelength of nuclei is smaller than typical interatomic distances. This is not the case 
for electrons which maintain the interatomic chemical bonds by the spatial extension of their quantum wavefunction. 
yf^ • Consequently, the motion of nuclei is essentially classical while electrons behave quantum mechanically. At low 
l/~s ' temperature, de Broglie's wavelengths are larger. For instance, in semiconducting nanodevices at subKelvin electronic 
^-Z ■ temperatures, the nuclei are nearly frozen while the electronic waves propagate ballistically in artificial circuits of a 
^— K I few hundred nanometers 0-3 ■ 

fvq r Newton's and Schrodinger's equations have in common to be time-reversal symmetric and, moreover, deterministic 

^— I in the sense that the time evolution they rule is uniquely determined by initial conditions belonging to either the 

;!^ , classical phase space of positions and velocities, or the Hilbert space of wavefunctions. In both schemes, the initial 
conditions themselves are not uniquely determined once for all, but free to take arbitrary values depending on the 
particular experiment considered. In this regard, the equations of motion describe all the possible histories that 
V^ ■ a system may follow, among which there exists the unique history actually followed by the system under specific 
conditions. Besides, most of the histories do not coincide with their time reversal in typical classical or quantum 
systems. For instance, the trajectory of a free particle with a non-vanishing velocity is physically distinct from its 
time reversal in the phase space of classical states: 

(rt,vt) = (vot-|-ro,Vo) ^ (r'^, v^) = (-Vot + Tq, -Vq) . (2) 

If the first trajectory was coming from the star Antares towards the constellation of Taurus on the opposite side of 
the celestial sphere, the second would travel in the opposite direction following a reversed history [5| . A similar result 
holds for the wavepacket of a free quantum particle. For such systems, the selection of initial conditions breaks the 
time-reversal symmetry. Indeed, the solutions of an equation may have a lower symmetry than the equation itself. 
This is well known in condensed matter physics as the phenomenon of spontaneous symmetry breaking. For example 
in ferromagnetism, the orientation of magnetization is determined by fluctuations at the initial stage of the process, 
i.e., by the initial conditions. 



These considerations extend to statistical mechanics where the time evolution concerns statistical ensembles de- 
scribed, in classical mechanics, by a probability density p evolving in phase space according to Liouville's equation: 

dtp^{H,p},, (3) 

where H is the Hamiltonian and {• , -jci the Poisson bracket or, in quantum mechanics, by a density operator p 
evolving according to von Neumann's equation: 

dtp^^[H,p] (4) 

where [ • , • ] is the commutator, i — %/— T, and fi. = /i/(27r) [6|. Both equations are symmetric under time reversal but 
they may admit solutions that correspond to histories differing from their time reversals. 

However, thermodynamic equilibria are described by Gibbsian statistical ensembles such as the canonical ensemble 
defined by the probability distribution 

p=lexp(-/3H) with /3-r^ (5) 

Z Kb -I 

for the temperature T and where fce is Boltzmann's constant. This distribution is a stationary solution of Liouville's 
or von Neumann's equation and it is time-reversal symmetric. As a consequence, the principle of detailed balancing 
holds according to which opposite fluctuations are equiprobable. The entropy production vanishes and there is no 
energy dissipation on average in systems at equilibrium. We notice that the thermodynamic equilibrium is a state 
which is stationary from a statistical viewpoint but dynamical from a mechanical viewpoint. 

During the last decades, remarkable advances have been achieved in the understanding of the statistical properties of 
nonequilibrium systems. These results find their origins in the study of large-deviation properties of chaotic dynamical 
systems, for which methods have been developed to characterize randomness in time (instead of space or phase space 
as done in equilibrium statistical mechanics) [3| . On the basis of dynamical systems theory, relationships have been 
established between dynamical large-deviation quantities and transport properties such as diffusion, viscosity, and 
electric or heat conductivities [8-18J. In this context, the further consideration of time reversal led to the discovery 
of different types of symmetry relations depending on the nonequilibrium regime. Indeed, systems may evolve out of 
equilibrium under different conditions: 

• In isolated systems, the solution of Liouville's equation is not stationary and evolves in time if the initial dis- 
tribution po differs from an equilibrium one. If the dynamics is mixing, the statistical distribution pt weakly 
converges towards a final equilibrium distribution after nonequilibrium transients. Transient states with ex- 
ponential decay exist in the forward or backward time evolutions, which are mapped onto each other by time 
reversal, while the equilibrium state is symmetric :14-16]. 

• In systems controlled by time-dependent external forces, the statistical distribution pt remains out of equilibrium 
since the relaxation towards equilibrium is not possible. For instance, the nonequilibrium work is a fluctuating 
quantity of interest in small systems such as single molecules subjected to the time-dependent forces of optical 
tweezers or atomic force microscopy jlft-.24:] . 

• In open systems in contact with several reservoirs at different temperatures or chemical potentials, a nonequi- 
librium steady state is reached after some relaxation time. Such states are described by stationary statistical 
distributions that are no longer symmetric under time reversal contrary to the equilibrium canonical state ([S]). 
Indeed, out of equilibrium, mean fluxes of energy or matter are flowing on average across the open system and 
this directionality breaks the time-reversal symmetry. In small systems or on small scales, the energy or particle 
currents are fluctuating quantities described by the nonequilibrium statistical distribution. 

Because the underlying microscopic dynamics is reversible, the fluctuations of nonequilibrium quantities obey 
remarkable relationships, which are valid not only close to equilibrium but also arbitrarily far from equilibrium [2l| - 
|35|. Several kinds of such relationships have been obtained for transitory or stationary nonequilibrium situations, for 
classical or quantum systems, and for Markovian or non-Markovian stochastic processes. 

Moreover, further time-reversal symmetry relations have also been obtained for the probabilities of the histories 
or paths followed by a system under stroboscopic observations |36l - l43j . It turns out that the time asymmetry of 
nonequilibrium statistical distributions implies that typical histories are more probable than their corresponding time 
reversal. In this regard, dynamical order manifests itself away from equilibrium (40[. Furthermore, the breaking of 
detailed balancing between forward and reversed histories is directly related to the thermodynamic entropy production, 
as established for classical, stochastic and quantum systems |36l - l43| . 



The purpose of the present chapter is to give an overview of the time-reversal symmetry relations established for 
nonequilibrium quantum systems |43|-|65|. Recently, important experimental and theoretical work has been devoted 
to transport in small quantum systems and, in particular, to electronic transport in quantum point contacts, quantum 
dots, and coherent quantum conductors [z-4|. In these systems, single-electron transfers can be observed experimen- 
tally and subjected to full counting statistics, allowing the experimental test of the time-reversal symmetry relations 
and their predictions. Here, the aim is to provide a comprehensive overview of the results obtained till now on this 
topic and to discuss some of the open issues. 

This chapter is organized as follows. In Section |lTl functional symmetry relations are established for quantum 
systems driven by time-dependent external forces in a constant magnetic field, which allows us to derive the Kubo 
formulae and the Casimir-Onsager reciprocity relations for the linear response properties. In Section lllli general 
methods are presented to describe open quantum systems in contact with several reservoirs and transitory fluctuation 
theorems are obtained. In Section lTVl an appropriate long-time limit is taken in order to reach a nonequilibrium steady 
state. In this limit, the stationary fluctuation theorem is obtained from the transitory one for all the currents flowing 
across the open system. In Section |Vl the current fluctuation theorem is used to obtain the response properties. The 
Casimir-Onsager reciprocity relations and fluctuation-dissipation relations are generalized from linear to nonlinear 
response properties. In Section IVTl these results are obtained for independent electrons in quantum point contacts or 
quantum dots. The generating function of full counting statistics is computed using Klich formula and the connection 
is established with the Levitov-Lesovik formula in the Landauer-Biittiker scattering approach. The equivalence with 
the Keldysh approach is discussed. In Section [VII[ the current fluctuation theorem is derived in the master-equation 
approach for the corresponding stochastic process. The second law of thermodynamics is deduced from the stationary 
fluctuation theorem for the currents and results on the statistics of histories are presented. In Section FVlIIi the general 
theory is applied to electronic transport in quantum dots, quantum point contacts, and coherent quantum conductors. 
Conclusions are drawn in Section [IXI 

II. FUNCTIONAL SYMMETRY RELATIONS AND RESPONSE THEORY 

We consider a quantum system described by the Hamiltonian operator H{t;B), which depends on the time t and 
the external magnetic field B. The time-reversal operator Q reverses the magnetic field but otherwise leaves invariant 
the Hamiltonian: 

QH{t;B)Q-^ ^H{t--B). (6) 

The reason is that the time-reversal symmetry of the electromagnetic interaction concerns the whole physical system 
including the external electric currents generating the magnetic field. The magnetic field changes its sign because the 
external currents are reversed under time reversal. At the level of the system itself, the symmetry should thus be 
implemented by the combination of the internal operator Q and the reversal of the background magnetic field. 

Since the system is driven by time-dependent forces^a comparison should be carried out between some protocol, 
called the forward protocol, and the reversed protocol |22h24| . 

In the forward protocol^ the system starts in the equilibrium distribution: 

-/3H(0;6) 

at the inverse temperature {3 = (fcET")"^ and the free energy F(0) = -kYiT\aZ{Q) with Z{Q) = tre~'^^(°'^). The 
system evolves from the initial time t — Q until the final time i = T under the unitary operator that is the solution 
of Schrodinger's equation: 

ih^UF{t-B)=H{t-B)Mt-B) (8) 

with the initial condition [/f(0;S) = 1. The average of an observable A is given by 

'iF(i))=tr/5(0)iF(t) (9) 

where 

AY{t)^ijl{t)AijY{t) (10) 



is the operator in the Heisenberg representation. 

In the reversed protocol, the system starts in the other equilibrium distribution: 

KT; -B) = ^^^^ (11) 

of free energy F{T) = —kBT\nZ{T) with Z{T) = tre"'^^^^'^®) and evolves under the reversed time evolution 
operator obeying 

ih^Unit; B) = H{T - t; B) Un{t; B) (12) 

at 

with the initial condition t/R(0; B) — 1. In the reversed protocol, the system thus follows the reversed external driving 
in the reversed magnetic field, starting at the time i = with the Hamiltonian H{T\ —B) and ending at the time 
t = T with the Hamiltonian H{0; —B). 
We have the 



Lemma: 51| The forward and reversed time evolution unitary operators are related to each other according to 

Unit;-B) = eUF{T-t;B)ul{T;B)e-^ with 0<t<T. (13) 

This lemma is proved by using the antiunitarity of the time-reversal operator, Eqs. (|S]) and p^ obeyed by the 
unitary evolution operators, as well as their initial conditions |5l| . 

Let us consider a time-independent observable A with a definite parity under time reversal, i.e., such that Q AQ~^ = 
ca-A with ca = ±1- The lemma allows us to relate its forward and reversed Heisenberg representations according to 

Apit) - eA llliT) e-i A^{T -t)Q Uy{T) . (14) 

In this setting, we have the 

Theorem: [51| If X{t) denotes an arbitrary function of time, the following functional relation holds: 

/q!o '** ^(*) ^f(*) Q-PH^i'T) ePHio) \ ^ e^'^'^-^ /e'-* ■1'^ '** ^C^"*) -^^^'^ \ (15) 

\ /f,B \ IK-B 

where Hp{7') ~ Up{'T) H{'T;B)tJF{'T) and AF — F{T) ~ -F'(O) is the free-energy difference between the initial 
eguilibrium states ill]) and f^ of the reversed and forward protocols. 

This theorem, which has been demonstrated in Refj5l|, extends results previously obtained in the restricted case 
where there is no change in free energy AF = [6a, [63 • 

Moreover, the quantum Jarzynski equality is deduced in the special case where A = in Eq. (ITSI) : 

In the left-hand side, the quantity inside the bracket can be interpreted in terms of the work performed on the system 
in the quantum scheme where von Neumann measurements of the energy are carried out at the initial and final times 
[44-46, 54-57]. In the classical limit where the operators commute, both exponential functions combine into the 
exponential of the classical work Wci = [-^f(T) — -ff (0)]^j and the classical Jarzynski equality [21| 

(e-^^"')F,B = e''^^ (17) 

is recovered for the nonequilibrium work performed on the system during the forward protocol. We notice that AF 
does not depend on the sign of the magnetic field because the canonical equilibrium distributions are time-reversal 
symmetric. 

Remarkably, the functional symmetry relation (J15p unifies the work relations and the response theory in a common 
framework. Indeed, the Kubo formulae as well as the Casimir-Onsager reciprocity relations can also be deduced 
from this relation. With this aim, we assume that the system is composed of N particles of electric charges {e„}^^i 
subjected to an external time-dependent electric field £,^{t) in the spatial direction v = 1, 2 or 3. The electric field is 



supposed to vanish at the initial and final times E^, (0) = E^, {T) = so that the free-energy difference is here equal to 
zero, AF — 0. The time-dependent Hamiltonian is given by 



N 



H{t) = %-£,{t)Y,enXn. (18) 

n=l 

with the position operators Xni> of the N particles. The observable we consider is the electric current density in the 
spatial direction /x = 1, 2 or 3: 

n— 1 

where V is the volume of the system. In order to obtain the linear response of the observable A with respect to 
the perturbation due to an electric field £v{t) of small amplitude, the functional derivative of Eq. ([T5|) is taken with 
respect to A(T) around A = 0, which yields 

A^[T)e-P"-^^U^"°) =/i\ . (20) 

/F,B \ /cq,B 

Developing this expression to first order in the electric field, the mean value of the current density at the final time 



t = T is obtained as [51 1 



3^{T) = Ia^{T)) = (i) + / dt £.iT~t) cl,^,{t;B) + 0{£l) (21) 



B \ ' cq.B Jq 

with the response function 

1 rP 



(j)^,{t-B)^- du{J,i-ihu)J^it))eq,B- (22) 

^ Jo 



Since the mean value of the current is vanishing at equilibrium (A) = 0, wc find the expression of Ohm's law 

\ / cq,B 

Jii = (^fj.iy £v (23) 

in the long-time limit T — > oo. The electric conductivity is given by 

/•CO 

<J^.u = / 0,..(i; B) dt = a^,{uj = 0; B) (24) 

Jo 



in terms of Kubo's formulae for the ac conductivities [68 1 



a^,{uj-B) = - dte'^* du{M~ihu)Jf,{t)),q,B- (25) 

^ Jo Jo 

These conductivities satisfy the Casimir-Onsager reciprocity relations [6^, [70] 

a^,{uj;B)^a,^{uj;~B). (26) 



In this way, the linear response properties are recovered from the functional symmetry relation (J15p . Higher-order 
terms in the expansion can be obtained as well for the nonlinear response properties. 

The functional time-reversal symmetry relation thus provides a unifying framework to study the fundamental 
properties of quantum systems driven out of equilibrium by time-dependent external fields. 

III. THE TRANSITORY CURRENT FLUCTUATION THEOREM 

A related problem of interest is the transport of particles and energy across an open quantum subsystem in contact 
with several reservoirs at different temperatures and chemical potentials. The subsystem is coupled to the reservoirs 
by a time-dependent potential V{t) during some time interval < t < T. The total system evolves in time under the 



quantum-mechanical unitary time evolution from the initial state specified by taking an equilibrium distribution in 
each separate part of the total system until the final time T, as schematically depicted in Fig. [T] Here, the general 
case is considered where the time-dependent potential may leave some permanent changes inside the subsystem and 
the reservoirs after the final time t = T. The flow of energy and particles between the reservoirs is determined by two 
successive von Neumann quantum measurements in each part of the total system, the first at the initial time i = 
and the second at the final time T. In this setup, a time-reversal symmetry relation is established for the transfers of 
energy and particles between the parts. 



Vit) 



m 



^ 



T 



FIG. 1: Schematic representation of tlie coupling of a subsystem to a pair of reservoirs by the potential V{t) during the time 
interval < t < T, for the protocol used to deduce the transitory current fluctuation theorem. At the initial time t = 0, the 
potential is vanishing, V'(O) = 0. At the final time t = T, the potential may be non vanishing, but then given by a sum of 
operators separately acting on the quantum state spaces of the subsystem and the reservoirs. 



The total Hamiltonian operator is given by 



H{t-B) = 



ns + EUn, 

^'s + ^3 = 1 '^j 



for t<0, 
V{t) for 0<i<r, 
for T <t, 



(27) 



where Hs and Hj denote the Hamiltonian operators of the isolated subsystem and the j^^ reservoir before the 
interaction is switched on. V(t) is the time-dependent interaction between the subsystem and the reservoirs. After 
the final time t — T, the subsystem and the reservoirs no longer interact with each other, but the final potential 
V{T) may have added possible contributions to the Hamiltonian operators of the subsystem and the reservoirs, Hg 
and H'j, which could thus differ from the initial ones. Moreover, the total Hamiltonian always obeys the time-reversal 
symmetry ([6]). 

The observables of the total system include not only the Hamiltonian operators but also the particle numbers of 
species a = 1,2, ..., c inside the subsystem, Nsa, and the reservoirs, Afja with j = 1, 2, ..., r. These latter observables 



are also time- reversal symmetric: QAfja O ^ = Sfja with j 
The total numbers of particles are given by 



S,l,2,... 



Na^Sfsa+^Sfja for a = l,2,...,c 
and they commute with the total Hamiltonian operator because they are conserved: 



[H{t;B),N^ 



0. 



(28) 



(29) 



The system may admit further constant numbers of particles if the reservoirs are not all connected to each other. 
If the total system was composed of d disconnected parts between which there is no flux of particles, the particle 
numbers of every species would be separately conserved in each part so that the system would admit c x d independent 
constants of motion. 

Moreover, before the initial time and after the final time, the reservoirs are decoupled from each other so that their 
Hamiltonian operator commutes with the particle numbers: 



'Hj,Nj'a 



for t < 



and 



h;,aG'c =0 



for r < t 



(30) 



for every j,j' — S, 1, 2, ...,r and a = 1, 2, ..., c. 

In nonequilibrium systems, the currents are determined by the differences of temperature and chemical potentials 
between the reservoirs. Since the reservoirs are large, the initial temperature and chemical potentials of the subsystem 



are not relevant if the subsystem is small enough. In this regard, we can simplify the formulation of the problem 
by regrouping the subsystem with one of the reservoirs, for instance the first one, and redefine the Hamiltonian and 
particle-number operators as follows: 

H,=Hs+Hi, H[=H's+K, Nio. ^ N'sa + Aa , for j = l and , . 

Hj = Hj , H'^ = H'j , Njc. = Mjc , for j = 2, ..., r ^ ' 



,r . 



In this case, the total Hamiltonian given by Eq. (1271) can be rewritten as 



Y.]=iHj for t<0, 

Ej=i Hj + V{t) for < t < r , (32) 



,E;=i^ for T<t, 

and the total particle numbers of species a=l,2,...,cas 






^a = >J^ja- (33) 



As in previous Section [TTl a comparison is made between the outcomes of some forward and reversed protocols. 
The forward protocol starts with the total system in the grand-canonical equilibrium state of the decoupled parts 
at the different inverse temperatures I3j — {k-B,T.j)^^ and chemical potentials /ijQ.: 

p(Q-B) = J|s,(6)-ie"^4^^(»)-5:=A...^..] = ■Qc-ft[^^(«)-^"'^-^^°-*^N (34) 

3 3 

where $j(S) = —k-QTj InSj(K) denotes the grand-canonical thermodynamic potential of the j**^ part in the initial 
equilibrium state. We notice that the grand-canonical potential is even in the magnetic field if the corresponding 
Hamiltonian has the time-reversal symmetry ([B]). The system evolves from the initial time i = until the final time 
t = T under the forward time evolution unitary operator obeying Eq. ([5]) with the initial condition Uf{0;B) — 1. 

In order to determine the fluxes of energy and particles, quantum measurements are carried out on the reservoirs 
before and after the unitary time evolution. 

An initial quantum measurement is performed which prepares the system in the eigenstate |^fe) of the energy and 
particle- number operators: 

t<0: 4|*fe) = £jfe|*fe), (35) 

Nja\-^k) = i^jafel^fc) . (36) 

After the time interval < t < T, a, final quantum measurement is performed in which the system is observed in 
the eigenstate |^J) of the final energy and particle-number operators: 

T<t: H'^\^[) = 41*;), (37) 

Nja\'i'[) = ^;„zl*;) . (38) 

We notice that semi-infinite time intervals are available to perform the initial and final quantum measurements of 
well-defined eigenvalues. Accordingly, this scheme based on two quantum measurements provides a systematic way 
to measure the energies and the numbers of particles transferred between the reservoirs during the time interval 
< t < T of their mutual interaction. During the forward protocol, the energy and the particle numbers are observed 
to vary in the j"^ part of the system by the amounts 

Ae, = e'^i - e,k , (39) 

l^Vjo, = V]^l - Vjo^k ■ (40) 

The probability distribution function to observe these variations is defined as 

PF[Ae,,Av,^-B) = Y,\{^[^^3-{^'3i-£3k)\ \{5[Av,^-{y'j^i-Uj^k)\ 

kl j ja 

xm{B)\UF{T;B)\^km\' {^kmrnB^km m 



in terms of Dirac distributions S{-). 

In the reversed protocol, the system starts in the other equilibrium distribution: 

3 3 

at the same inverse temperatures /3j = {k^Tj)~^ and chemical potentials /Xj-q as in the forward protocol. Here, 
^'j{B) — —ksTj lnS^(;B) denotes the grand-canonical thermodynamic potential of the j*^ part in the final equilibrium 
state and reversed magnetic field. The reversed time evolution unitary operator obeys Eq. (I12p and it is related to 
the one of the forward protocol by Eq. p^ with t — T oi the Lemma in Section |lll Before and after the reversed 
protocol, quantum measurements are similarly carried out in order to determine the variations of energy and particle 
numbers in the reservoirs. The probability distribution function of these variations during the reversed protocol is 
defined by an expression similar to Eq. (|4T1) . 

Comparing the probability distribution functions of opposite variations during the forward and reversed protocols, 
the following symmetry relation is obtained 



PR(-Ae„-Aiyj„;-S) 



with the differences A$j = <I>^- — ^j in the grand potentials |52l.[58|. 

For the forward and reversed protocols, the generating functions of the statistical moments are defined as 

Gf,r{^j,V3o.;13) = / n ^^^J ^^^J" e" ^' «^^^^-^° "-^''^ Pf3(A£,, Ai.,„;B) (44) 

ja 

in terms of the so-called counting parameters {£,j,'rija}- The statistical moments of the energy and particle-number 
variations can be obtained from this generating function by taking derivatives with respect to these counting param- 
eters. Now, the symmetry relation (j43l) can be equivalently expressed as 

gF{^j,V3c.;B) = e- ^. ft^*^ GniPj - (j, -Pjt^ja - V3c.;-B) (45) 

in terms of the temperatures and chemical potentials of the reservoirs. We point out that this relation is only 
symmetric with respect to the inverse temperatures and the chemical potentials of the reservoirs. 

Such symmetry relations are useful for quantum system driven by time-dependent external forces. If the relation 
(I43p is restricted to the energy variation, the fluctuating quantity is the work W performed on the system and we 
recover the quantum version of Crooks' fluctuation theorem [23| 

Pf{W;B) _,^w-af) 



Pn{-W;-B) 



^mw-i^n (4g) 



with the free-energy difference AF ^ F' — F. We notice that this relation implies Eq. p^ [56]. Similar symmetry 
relations may be obtained for cold atoms or molecules in rotating frames, in which case the rotation rate fl plays the 
role of the magnetic field B. 

IV. FROM THE TRANSITORY TO THE STATIONARY CURRENT FLUCTUATION THEOREM 

If r infinitely large reservoirs are coupled together via the subsystem of interest, a stationary state will be reached 
in the long-time limit |7ll473l |. According to thermodynamics, the whole system is at equilibrium if the temperature 
as well as the chemical potentials of the different particle species are all uniform, i.e., if every reservoir shares the 
same temperature and chemical potentials. This is no longer the case if the reservoirs have different temperatures 
or chemical potentials whereupon energy and particles are exchanged between the reservoirs across the subsystem 
coupling them together. Therefore, energy and particle currents are induced by the so-called thermodynamic forces 
or affinities defined in terms of the differences of temperatures and chemical potentials with respect to some reference 
values [74.-76]. 

If we consider an open subsystem in contact with r = 2 reservoirs as depicted in Fig. [51 the thermodynamic forces 
or affinities are defined as 

Ae = (3i-(32, (47) 

Aa = |32^i2a - |3l^J-la for a = I,2,...,c, (48) 



N 



I — I Z3 



V(t) 



U 



T-t, 1 



FIG. 2: Schematic representation of the coupling of a subsystem to a pair of reservoirs by the potential V{t) during the time 
interval < f < T, for the protocol used to obtain the stationary current fluctuation theorem. Here, the time-dependent 
potential has the symmetry V{t) — V{T — t), remaining constant during the time interval to < t < T ~ to, and vanishing at 
the initial and final times: V'(O) — V{T) = 0. 



driving respectively the energy and particle currents from reservoir No. 2 to reservoir No. 1. The nonequilibrium 
conditions are thus specified by these affinities. 

At the microscopic level of description, this system is described by the following Hamiltonian operator: 



H{t; B) ^ Hi + H2 + V{t) for < t < T 



(49) 



Here, the interaction V{t) is switched on during a short time interval < t < to at the beginning of the forward 
protocol with to <C T, remains constant V{t) = Vo during most of the protocol ior to < t <T—to, and is switched off 
at the end when T — to < t < T , a,s shown in Fig. [51 Moreover, the time dependence of the interaction is assumed to 
be symmetric under the transformation t ^ T —t. Therefore, the forward and reversed protocols are identical except 
for the reversal of the magnetic field. Consequently, A^j = and the functions ([H]) of the forward and reversed 
protocols coincide, Qp = ^r = Q-j-. In this case, the symmetry relation (j45p becomes 



GT{£,l,^2,Vla,V2a;B) = Gril^l -S.l,h -6,-/3lAila -?7la,-^2M2a - V2a; -B) 



(50) 



For applying to nonequilibrium steady states, the shortcoming is that this relation is not a symmetry with respect to 
the affinities or thermodynamic forces driving the system out of equilibrium. Accordingly, this relation does not yet 
concern nonequilibrium steady states associated with given affinities and further considerations are required. 
We introduce the cumulant generating function as 



lim 



T 



In Qr ■ 



According to Eq. (|44p . the statistical averages of the energy and particle-number variations are given by 

(Ae,) , dQ 



dQ 



lim 

7~— -^00 



T 



and 



drija 



= lim i^ 

r^oo T 



(51) 



(52) 



while further differentiations give the cumulants. Besides, the short time interval to is assumed to be constant in the 
long-time limit T — > 00. Remarkably, it is possible to prove that 

Proposition: [52l | // the long-time limit i51]) exists, the cumulant generating function only depends on the differ- 
ences C = ^1 ^ C2 cind rja — rjia — ?/2a between the counting parameters associated with the reservoirs: 



Q{S„ria;B) = ^lini^ -— In Qr Uo + ^'^o - i^^lQa + y'^^Oa - -^'^^ 

and the cumulant generating function of the energy and particle currents is given by 

Q{i,Va;13)^ lim -^lngt{^,7^^■,B) 
t-^00 t 



for any {^o,%a} 



(53) 



(54) 
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where 

Gti^, r,^;B) = tr po e'"' e'^^"^-^'^^-^^ ^(*i=-^2.) ^-^m g|(//i-ff.)+E„ ^(A.-n,^) (55) 

with t = T — 2to, H = Hi + H2 + Vq, and po is the initial grand- canonical distribution ^34^ , which fixes the 
temperatures and the chemical potentials of the reservoirs. 

This proposition is proved by bounding together the functions (|50p and ([55|l according to 

Lgt<gr<Kgt (56) 

where Qt — Gt (Ci — C2 : 'yia — rj2a ; B) and the two factors L and X arc independent of the time interval 7" [52| . Therefore, 
the contributions of the factors L and K disappear in the long-time hmit T ~> 00. Moreover, hm7-_j.oo (V^) = 1 since 
t = T — 2to and to is constant. For these reasons, the resuUs (|53l) and (l54l) are obtained. Q. E. D. 

The inequahties (|56p are established in several steps pz] . The first two steps consist in removing the contributions 
from the initial and final lapses of time to, during which the interaction is switched on or off. Indeed, these contributions 
become negligible in the long-time limit. The third step has the effect of letting appear the differences ^1 — ^2 
and Tjia — ri2a between the counting parameters of the two reservoirs, which is the main result of the proposition. 
Accordingly, Eqs. ((M)) and ([55)) genuinely concern the fluxes of energy and particles transferred between the reservoirs. 

Now, the symmetry relation (j50p implies the following 

Current fluctuation theorem: 52] As the consequence of the time-reversal symmetry Q H{B)Q^^ = H{—B), 
the cumulant generating function \54^ of the energy and particle currents satisfies the symmetry relation: 



Q{LVc.\B)^Q{Ae-£.,A^-7^^--B) (57) 

with respect to the thermal and chemical affinities \41^ and C^. 

We notice that the function (I55|) takes the unit value if ■^ = ?7q = so that the cumulant generating function 
vanishes with the counting parameters: 

Q{Q,Q-B)^Q{Ae,A^--B) = Q. (58) 

The current fluctuation theorem (|57|) can be extended to the case of systems with more than two reservoirs [52] . 
If all the r reservoirs are coupled together, the total energy as well as the total numbers of particles are conserved 
in the transport process so that there are {c -\- 1) constants of motion. The grand-canonical equilibrium state of 
each reservoir is specified by one temperature and c chemical potentials. One of the r reservoirs can be taken as a 
reference with respect to which the affinities are defined. Consequently, the nonequilibrium steady states are specified 
by P = (c -f l)(r — 1) different affinities and so many independent currents may flow across the open subsystem in 
between the reservoirs. 

Some systems may be composed of several separated circuits between which no current is flowing^ ^This is the 

case for instance in quantum dots monitored by a secondary circuit with a quantum point contact |59l - l61| . Both 
circuits are only coupled by the Coulomb interaction so that there is no electron transfer between them. In such 
systems, there are more independent affinities and currents. If a system composed of r reservoirs is partitioned into 
d disconnected circuits each containing at least two reservoirs, a reference reservoir should be taken in each separate 
circuit. Therefore, (r — d) reservoirs are controlling the (c -t- 1) possible currents. The nonequilibrium steady states 
are thus specified hy p — (c+ l)(r — d) different affinities and so many independent currents can flow in such systems. 

V. THE CURRENT FLUCTUATION THEOREM AND RESPONSE THEORY 

If we collect together all the independent affinities and counting parameters as 

A^{AjE,Aj^} (59) 

A = (0 - S.k,Tllja - Vka] (60) 

where k stands for the indices of the d reference reservoirs, the current fluctuation theorem can be expressed as 

S(A,A;S) = Q(A-A,A;-6) (61) 



11 

where we have introduced explicitly the dependence of the generating function on the affinities A which specify 
the nonequilibrium steady state of the open system. The maximum number of independent affinities and counting 
parameters is equal to p — {c + l){r — d) . If the total system is isothermal, there is no energy current and the number 
of independent affinities and counting parameters is equal to p = ex (r — d) where c is the number of different particle 
species. 

All the statistical cumulants of the energy and particles transferred between the reservoirs are obtained by taking 
successive derivatives with respect to the counting parameters and setting afterwards all these parameters equal to 
zero. The mean values of the currents are given by the first derivatives, the diffusivities or second cumulants by the 
second derivatives, and similarly for the third or higher cumulants: 

MA;B) = |^(0,A;S) (62) 



with a,l3,... — 1,2,..., p. All these cumulants characterize the full counting statistics of the coupled fiuctuating 
currents in the nonequilibrium steady state associated with the affinities A in the external magnetic field B. 

At equilibrium, the mean currents vanish with the affinities. Therefore, we may expect that the mean currents can 
be expanded in powers of the affinities close to equilibrium: 



J„ == ^ L^^p Ap + -Y, Ma,P-f ^/3 A + • • • (65) 

which defines the response coefficients: 



2 

;3 /3,7 






We observe that the cumulants are given by successive derivatives with respect to the counting parameters. On the 
other hand, the response coefficients are given by one derivative with respect to a counting parameter and further 
successive derivatives with respect to the affinities. The remarkable result is that differentiating the symmetry relation 
(j6ip with respect to the affinities leads to derivatives with respect to the counting parameters. In this way, the current 
fluctuation theorem implies fundamental relationships between the cumulants and the response coefficients. 

Considering all the second derivatives of the generating function with respect to the counting parameters and the 
affinities, it is possible to deduce the Casimir-Onsager reciprocity relations 

L„,^(6) = Lp^^i-B) , (68) 

as well as the identities 

L^Al3) + Lp,^iB) = 2D^piO;B) (69) 

which relate the linear response coefficients to the diffusivities at equilibrium where A = 0. 
In the absence of magnetic field S = 0, we recover the Onsager reciprocity relations 

ia,/3(0)=i/^,„(0), (70) 

as well as the expression of the fluctuation-dissipation theorem according to which the linear response coefficients are 
equal to the diffusivities of the currents around equilibrium: 

L^,p{0) ^ Dafi{0;0) . (71) 



12 



Going to the third derivatives of the symmetry relation (|6ip . the third cumulants (|64)) characterize a magnetic 
asymmetry at equihbrium where A = and they are related to the responses of the diffusivities with respect to the 
affinities according to 



C„^^(0;6) 



-a,,(0; -S) = 2 ^(0; S) - 2 ^(0; -S) . 



dA^ 



dA^ 



(72) 



Moreover, the following relations hold for the nonlinear response coefficients at second order in the perturbations with 
respect to equilibrium [52 1 



M„,^^(6) + M„,0^(-S) 



2^^(0;6) + 2 



dA. 



^rn.r,, , o^^"7(o._g)^ 



M^j^iB) + Mii,^^{B) + M^,„(i{B) = 2 



dAp 

/dDj3y dDja dDai3 



V dAc, dAf3 



dA^ 



2 C'a,37 



(73) 
(74) 



A=0:B 



In the absence of magnetic field B = 0, the magnetic asymmetry disappears because the third cumulants vanish, 
Cap-yiO', 0) — 0, and the response coefficients are given in terms of the diffusivities as follows |77h79| 



M„,0^(O) = 



I dA^ 



dAt 



P / A=O;e=0 



(75) 



These relations are the generalizations of Eq. (I7T|) to nonlinear response coefficients. 

Such relations, which can be extended to higher orders as well [52l l77H79l |. find their origin in the microreversibility 
expressed by the current fluctuation theorem (|6ip . 



VI. THE CASE OF INDEPENDENT PARTICLES 



In order to illustrate the previous results, let us consider systems with independent fermionic particles ruled by the 
Hamiltonian: 



H 






(76) 



written in terms of the the anticommuting fleld operators ■(/'^(r) and the one-particle Hamiltonian operator: 

fi2 



-^V^+.(r) 



with 



r = {x,y,z) . 



(77) 



The confining potential u(r) is minimum in the conducting region of the circuit and presents a barrier at the borders 
of this region. Such confining potentials can be fabricated by nanolithography on semiconducting heterojunctions and 
various shapes can be given to the circuit 3]. If the temperature is of the order of a few Kelvins or less, only the 
electrons around the Fermi energy are transported. The de Broglie wavelength can be larger than the interwall 
distance so that the circuit behaves as an electronic waveguide. 



(a) 



(b) 




FIG. 3: Schematic representation of different electronic circuits: (a) quantum point contact between two reservoirs; (b) quantum 
dot in tunneling contact with two reservoirs. 



A quantum point contact (QPC) is a circuit with a barrier separating two semi-infinite waveguides, which form 
two reservoirs (see Fig. |3^). The barrier is a bottleneck for the transport of electrons from one reservoir to the other. 
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This bottleneck presents a saddle point in the potential energy landscape. If the Fermi energy is lower than the 
energy of the saddle point, the transport proceeds by quantum tunneling through the barrier. If the Fermi energy is 
larger, the electronic waves undergo direct scattering on the obstacle formed by the barrier. Different models can be 
envisaged for the potential of the QPC. At large distance from the barrier, electrons can be assumed to propagate in 
a waveguide with a potential u(r) = uj_{y,z), which is invariant under translations in the a;-direction and confining 
in the transverse y- and z-dircctions. The propagation modes in the infinite waveguide are given by the eigenstates 
of the one-electron Hamiltonian: 

ho i'k^nyn, (r) = £k^nyn, V'fc.n^n, (r) (78) 

where k^ is the wavenumber in the x-direction of propagation while Uy and n^ are the quantum numbers of the modes 
in the transverse directions. Each of these modes is a possible channel of conduction, which opens at the energy 
threshold given by the minimum of the energy band Sk^n n^ at kx — 0. 

The propagation modes ^pk^n n,(r) are scattered by the barrier. If the potential can be assumed to be the sum of 
the transverse and barrier potentials, u(r) — u^\{x) + u±{y,z), the problem is unidimensional for the passage of the 
barrier. If the incoming wavefunction is a plane wave in the positive x-direction, the scattering generates reflected 
and transmitted waves: 

f ^= (e*''^ + Tfc e"*''''^) for x <0 , 

^"^ ^ 1 ^ij.e**-^ for X > , ^ ' 

where r^ and t^ denote the reflection and transmission amplitudes while fc = fc^, is the wavenumber. At the energy 
e = £{k), the outgoing waves are related to the incoming waves by the unitary scattering matrix according to 

*out = S{e) *in (80) 

where Vl^in and vE'out denote vectors with two components. The first component is the wave amplitude in the left-hand 
reservoir and the second in the right-hand reservoir. Therefore, the unitary scattering matrix of the one-dimensional 
barrier is given by 

S{e) = ( \^ ^^ \ such that S\e) S{e) ^ 1 . (81) 

\^k '^k J 

The unitarity of the scattering matrix expresses the conservation of probability and, in particular, the fact that the 
transmission probability T(e) = |ifcp and the refiection probability R{e) = jr^p add to the unit value: r(e)-f i?(e) = 1. 
Here, the scattering matrix is 2 x 2 for every transverse mode because the potential defines a problem that is separable 
into the longitudinal wave and the transverse modes. This would not be the case if the potential was no longer the 
sum of two potentials in perpendicular directions. In more complicated circuits such as billiards, it is known that the 
scattering matrix is infinite and couples together the different transverse modes [80|, ISlI • 

In the case of separable potentials, every mode is scattered independently of the other ones and the scattering is thus 
characterized by the transmission probability T{e) — |ifcp, which depends on the wavenumber k — k^ or, equivalently, 
on the corresponding energy s = e{k) of the electron in the longitudinal direction. A well-known example is the 
inverted parabolic barrier u(x) = uq — mX^x^/2, for which the transmission probability is given by j4!. i82i| 



T{s) 



1 + exp -27r — — — 



(82) 



The transmission probability converges to the unit value at high energy well above the height of the barrier. For 
energies lower than the height of the barrier, the transmission proceeds by tunneling. The broader the barrier, the 
smaller the transmission probability. Since the inverted parabolic potential is unbounded from below, the transmission 
probability (j82p remains positive for e < and vanishes only in the limit e — > — oo. For potentials which are vanishing 
at large distances, the transmission probability is strictly equal to zero if e < 0. 

A further example of electronic circuit is composed of a quantum dot between two reservoirs, as depicted in Fig.[3}3. 
A quantum dot is an artificial atom in which several electrons are confined on quasi discrete energy levels [3, |j] . In 
general, there exist several electronic energy levels but the transmission can be observed around a single electronic 
level. If the tunneling to the reservoirs is small enough, the energy levels are observed as resonances in the transmission 
probability. For quantum circuits such as classically chaotic billiards, the scattering resonances may form an irregular 
spectrum in the complex plane of the wavenumber or the energy and semiclassical methods become useful [80|, |8l| . 

An important remark is that the Hamiltonian operator (1761) can also be written in the form: 

H^Hi^ + Hn + Vo (83) 
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where Hi^ and H^^ are the Hamiltonian operators of the left- and right-hand reservoirs, while Vq is the interaction 
between the reservoirs, this interaction being at the origin of the scattering. 

For these systems, we consider the function (|55|) in the absence of magnetic field and for one species of particles, 
namely, the electrons. This function has the following form: 



with 



and 



^t(C,r?)=trpoe-'^'e^° 



i=|(7?R-7^L) + |(iVR,-7VL) 



Af = e'"' Ae-'"' 



(84) 



(85) 



(86) 



where Hi = Hr and H2 = i?L, in order to define the affinities (j47|) and (|48|) corresponding to currents from the 
left-hand reservoir to the right-hand one. 

According to Eq. (|M)). the initial density operator is given by 



-B 



PO 



with 



tre" 



B = /3L(-ffL - MLiVL) -t- l3n{Hn - m^n) 



(87) 



in terms of the inverse temperatures, the chemical potentials, and the particle-number operators, TVl and TVr, for 
both reservoirs. Since [A, B] = 0, the function (|84p can be written as 



Gt{^,v) 



tre" 



-At „A-B 



tre 



In systems with independent particles, many-particle operators such as the Hamiltonian and particle-number op- 
erators are of the form 



X = ^Xkiclci = T{x) 



(89) 



kl 



where x — {xki) denotes the corresponding one-particle operator, while cj, and Ck are the anticommuting creation- 
annihilation operators in the one-particle state k. The correspondence between one-particle and many-particle oper- 
ators is thus established by the mapping F defined by Eq. (|89|). In this framework, Klich has shown that the trace of 
products of exponential functions of many-particle operators can be expressed as appropriate determinants involving 
the one-particle operators [8J, |8J|. For fermions, Klich's formula reads 



tr e^ e^ = tr e^^*) e^^^^ = det (l + e* e^) 
Applying Klich's formula to the function ([55]) . this latter becomes 

-1 



gti^,ri)=dct 



1 + c 



-b 



1 + c-'^' e"-'' 



= det 



1 + / (e-°* e° - 1 



(90) 



(91) 



with the operator: 



/ 



1 



3;9L(hL-/JL"L)-|-/5R('»R-MR"R) -|- 1 



of the Fermi-Dirac distributions in the reservoirs. 

In the long-time limit, the unitary scattering operator is defined as 

S = lim e''">*/^ e-'^* e'^'>*^^ 



(92) 



(93) 



in terms of the full Hamiltonian h and the non-interacting Hamiltonian Iiq = /il + ^r which commutes with the 
operators a and b. Since 



dt = e'''*ae-'''' 



(94) 
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we have that 

in the long-time hmit. Using the cominutativity of ho with a and /, the function (|9ip can be written as 

^t (e, ?7) = det [l + / (^t e-a 5 e^ - l) ] (96) 

For this formula to make sense, we have to notice that the determinant should be taken over an appropriate discrete 
set of electronic states forming a quasi continuum in the long-time limit. Indeed, the quantity in the left-hand side 
of Eq. ([M)) concerns the whole many-particle system although the expression in the right-hand side concerns single 
independent electrons flowing one by one across the system. Every one-particle operator can be decomposed on the 
eigenstates of the non-interacting Hamiltonian /iq = ^l + ^r- This is the case in particular for the scattering operator: 

deS{e)d{e~ho) (97) 

where S{e) is the 2x2 scattering matrix (|8ip acting on the two wave amplitudes with opposite wavenumbers ±fc 
at the given energy e = e{k). Accordingly, the function (P5)) becomes a product over all the relevant single-electron 
states {£,cr} of the determinants of the corresponding operators at this particular energy: 

Gtii, v)=Y[dcti^l + /(e) [^t (e) e-^^^) ^(e) e^^^^ - l] } (98) 

£,fT 

where the matrix containing the Fermi-Dirac distributions of the left- and right-hand reservoirs is given by 

/(e) = ( \(^) \] with /l r(£) = , , ^ — ^ (99) 

and 

In order to define the cumulant generating function with Eq. (IM)) . the logarithm is taken of the function (P5|. The 
logarithm of the product over the relevant states {e, ct} gives a sum of logarithms. For a process lasting over the lapse 
of time i, the spacing between the relevant energies {e} is equal to Ae = 2Trh/t. Since these relevant energies form a 
quasi continuum, the sum is replaced by an integral in the long-time limit. Moreover, the spin orientation takes the 
two values ct = ± for electrons. For spin one-half particles such as electrons, we thus find that 



t^oo t ^^' J 2T:h^' 



(101) 



where 2s denotes the factor two due to the electron spin Q. 

After evaluating the determinant, the cumulant generating function (1541) is finally given by 

Q(C, ^7) = -2s / ^ In {1 + T{e) [/l(1 - /r) (e-^^-" - l) + /r(1 - /l) (0^?+" - l)] } (102) 

for independent electrons. Therefore, we have recovered the Levitov-Lesovik formula for the full counting statistics of 
electron quantum transport |85l l86l|. The symmetry of the fluctuation theorem is satisfied: 



Q(C,^)-Q(A£-C,^iv-'7) (103) 

with respect to the thermal and chemical affinities 

^B = /3r - /5l , (104) 

An = /3l ml - /3r MR , (105) 

defined by Eqs. (|Tf|) and (|^51) if the reservoir No. 1 (respectively No. 2) is identified with the right-hand (respectively 
left-hand) reservoir. 
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FIG. 4: The cumulant generating function (|102|) of the particle current versus the counting parameter r] for the inverted 



parabolic barrier u{x) 



Mo 



mX X /2 with the transmission probability 



The temperature difference and, thus, the 



thermal affinity (|104p are vanishing. In this regard, the generating function is considered at ^ = 0. The chemical potentials of 
the left- and right-hand reservoirs are given by /il.r = /io ± eV/2 with eV = 1 and po = 0, 0.5, 1, 1.5, 2. The other parameter 
values are uq = 1, h\ — 3, and /3 = {ksT)^^ — 1. Planck's constant is taken as h = 2s. The plot confirms the symmetry 
•q —7- yljv — 77 of the fluctuation theorem with respect to the chemical affinity Am = /3(ml — Mr) = lieV = 1. 



The cumulant generating function (|102p with ^ — Ae = is depicted in Fig. |4] for the inverted parabolic barrier 
with the transmission probability ([82]). We see that the generating function has indeed the symmetry 77 — > An — 77 of 
Eq. (|103l) . The generating function is depicted for several values of the mean chemical potential //q = (ml + Mr)/2 of 
both reservoirs. We observe that the generating function increases with the mean chemical potential /xg because the 
transmission probability (|82|) does so. 

The generating function given b y th e Levitov-Lesovik formula can be interpreted in terms of an exclusion process 
ruled by a trinomial distribution [83]. Suppose that, at every relevant energy, three types of random events may 
happen for an electron: it may move from the left- to the right-hand reservoir with the probability p+, in the other 
direction with the probability p- , and stay in some reservoir with the probability p^ — 1 ~ p+ — p_ . The probability 
that A'o electrons do not move, N+ electrons move to the right-hand reservoir, and N- to the left-hand reservoir 
during the time interval t is given by the trinomial distribution: 



P, 



N\ 



NoN+N- 



No N+ N- 
Po P+ P- 



NolN+lN^l 

For this random process, the generating function is defined as 

1 



with 



N = No + N+ + N^ 



Qiv) 



lim 



t 



ln(e-"^^) 



(106) 



(107) 



where AA^ = A^+ — N^ is the total number of electrons that have been transported from the left- to the right-hand 
reservoir during the time interval t. Taking the average over the trinomial probability distribution, we find that 



Q{r]) = -r In [l + p+ (e-" - 1) + P- (e" - 1)] 



(108) 



where r = \mit^c>o{N/t) is the attempt frequency, i.e., the mean number of attempted electron transfers per unit time 
[j] . We observe that the special form of the generating function (|102|) is indeed obtained at every relevant energy if 
we take p+ = T{e)fi^{l — /r) and p_ = T(e)/R(l — /l), as it should. The generating function of the trinomial process 
has the symmetry 



Qiv) ^ Qi^N - v) with An = In 



P+ 
P- 



(109) 



This provides a stochastic interpretation for the quantum transport process described by the Levitov-Lesovik for- 
mula pn^ . 
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Many known results for quantum transport can be deduced in this framework. In particular, the mean currents are 
given by 

Jn = ||(0, 0) = 2, 1 ^ T{e) (/l - /r) , (110) 



Je = ^{0,0} ^ 2, 1 ^eT{s){f^-M) 



(111) 



which are the well-known Landauer formula for the particle current and the equivalent formula for the energy current 
[88j . On the other hand, the noise is characterized by the diffusivities: 

with 

i?™ = y'^£"n£)[/L(l-/R) + /R(l-/L)-T(e)(/L-/R)2] for m = 0,l,2 (113) 

as reported in Refs. 0, |8^, |90|- The power spectrum of the electric noise is defined as 

S{u;) = / dte"^'{Ai{t) A/(0) + A/(0) A/(t)) (114) 

with AI{t) — I{t) — {I{t)) 90]. Since / — gJat, the zero-frequency limit of the noise power is thus related to the 
corresponding diffusivity by 

S{uj = 0) = ie'^DNN ■ (115) 

Close to equilibrium, the mean currents can be expanded in powers of the affinities (|104l) and (|105p . However, 
alternative expansions can be considered in terms of the potential difference V and the temperature difference AT 
between the left- and right-hand reservoirs at the temperatures Tl,r = Tq ± AT/2 and chemical potentials /iL,R — 
fio ± eV/2. If we introduce the electric and heat currents as 

I^eJN, (116) 

Jq = Je - ^J.oJN , (117) 

with the electric charge e, the alternative expansions around equilibrium read: 

AT 
I^GV + Li-^, (118) 

-to 

AT 
Jq^LiV + L2-^, (119) 

-to 

if higher-order terms are neglected. In this linear approximation, G = Tq is the electric conductance, Ti the thermo- 
electric coupling coefficient, and L2 the thermal conductance, which are given by [88| 



L^^2,e'-'^ J ^ie^fiorT{s)f~^) with m = 0,l,2 



(120) 



We notice that the Onsager reciprocity relation is satisfied. 

As the example of Eq. (j82p shows, the transmission probability tends to the unit value for energies higher than 
the barrier in every open channel corresponding to a propagation mode. In this limit where T(e) — 1, the electric 
conductance given by Eq. (|120l) with m = takes the universal value: 



2^_ 2e_2 
2T:h ^ "/T 



G. = ^--^- (121) 



If several channels are open at the Fermi energy ep , the zero-temperature conductance is given by 

2e 
1 



^=¥E^«(^f) (122) 
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where T„(£) is the transmission probabiUty of the channel corresponding to the mode n = {ny,nz) in Eq. ((78|) . As 
the Fermi energy increases, the channels open successively and their transmission probability tends to the unit value, 
which explains the phenomenon of conductance quantization [2-0, [20, ISlI ■ 

In summary, the results obtained in the previous sections allow us to deduce systematically many well-known results 
on quantum transport such as the Levitov-Lesovik and the Landauer formulae in the scattering approach, which is 
appropriate to treat coherent quantum transport. 

Similar considerations also apply to boson transport [49l l53j as well as to many other quantum transport processes 
in physics or chemistry where the theory of full counting statistics extends and complements standard scattering 
theory and reaction rate theory. 

We notice that these quantum transpo rt p roperties can also be deduced in the so-called Keldysh formalism, which 
has several conceptual advantages [3|, Ij, |92| |. On the one hand, it formulates transport in terms of quantum fields, 
allowing in principle to treat the spatial dependence of local properties such as densities or current densities. On 
the other hand, the Keldysh formalism provides the systematic introduction of all the moments of the quantum 
fields, starting from the nonequilibrium Green's functions, which is suitable for perturbative calculations. Several 
nonequilibrium Green's functions are introduced corresponding to different time contours, which come with the 
doubling of the space of quantum states from pure states to statistical mixtures described by density operators. An 
equivalent theory is the Liouville-space formalism developed by Mukamel and coworkers [58l l93l |94| , which starts from 
the quantum Liouville superoperator, i.e., the generator of von Neumann's equation Q. Another equivalent theory 
is provided by the formalism of thcrmofield dynamics, which also takes into account the need to double of the state 



electron transport in molecular junctions |96l. |97| 



space for the description of nonequilibrium systems |95| . The latter formalism has recently been applied to study 



VII. TIME-REVERSAL SYMMETRY RELATIONS IN THE MASTER-EQUATION APPROACH 

A. Current fluctuation theorem for stochastic processes 

The master-equation approach allows us to establish connections with the theory of stochastic processes. The first 
step of this approach consists in identifying the relevant coarse-grained states of the system or subsystem of interest. 
The second step usually proceeds with a perturbative calculation with respect to some small coupling parameter A, 
which typically controls the time scale separation between the slow time evolution of the coarse-grained states and 
the fast dynamics of the other degrees of freedom. The Hamiltonian may have the form 

H = Ho + XV (123) 

where the unperturbed Hamiltonian Hq leaves invariant the coarse-grained states and the perturbation XV causes the 
interaction between the coarse-grained states and the other degrees of freedom. Without a clear separation of time 
scales, the time evolution of the coarse-grained states continues to keep the memory of the past and is non-Markovian. 
In quantum-mechanical systems, non-Markovian behavior can manifest itself as a slippage of initial conditions on the 
time scale of the fast degrees of freedom in the early time evolution of a subsystem in contact with a reservoir. The 
slippage of initial conditions allows the density operator to remain positive definite [93, ll^]- Thereafter, the density 
operator of the subsystem follows a time evolution which is essentially Markovian on the long time scale characterizing 
the interaction of the coarse-grained states with the fast degrees of freedom. We notice that, although the subsystem 
density operator evolves slowly, the individual realizations present quantum jumps between the coarse-grained states. 
These jumps occur over the short time scale of the fast degrees of freedom and the dwell times between the jumps are 
of the order of the long time scale. Indeed, these dwell times are inversely proportional to the transition rates, which 
are of the order of A^ according to second-order perturbation theory. 

In transport problems where the subsystem of interest is coupled to several reservoirs, it is important to specify 
the coarse-grained states not only of the subsystem but also of the reservoirs, in order to proceed with the count- 
ing of particles and energy transferred between the reservoirs. This corresponds to taking two successive quantum 
measurements, as done in Section IlIII and this is sometimes referred to as unraveling the master equation [100l4l02 | 



Consequently, the master equation rules the density operator of the subsystem conditioned to the numbers of particles 
and energy that have been transferred between the reservoirs since the initial time. 

Many experimental observations are not sensitive to the quantum coherences described by the off-diagonal elements 
of the density operator. Under such circumstances, the knowledge of the diagonal elements, i.e., the probabilities of 
the coarse-grained states are enough for the description of the observations. In these cases, the time evolution of the 
probabilities is ruled by a Markovian mas ter e quation. For these stochastic processes, a current fluctuation theorem 
has been proved using graph analysis [78l , Il03l | . This theorem is expressed by the symmetry relation (I6ip in terms of 
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the cumulant generating function of the fluctuating currents. Another proof of this theorem is based on the symmetry 
of the master equation modified to include the counting parameters, as explained here below. 

We consider electron transport through an isothermal open subsystem in contact with r reservoirs. If the circuit is 
composed of d disconnected subcircuits, the numbers of electrons in these subcircuits are so many conserved quantities. 
We suppose that each subcircuit connects at least two reservoirs so that r > 2d. In this system, the nonequilibrium 
conditions are imposed with p = r — d differences of chemical potentials, which define the independent affinities 
Aj = P{fij — fik) with j = l,2,...,p and k = l,2,...,(i. These affinities may drive so many independent currents. 
We denote by n = {nj}'^^^ the random numbers of electrons transferred from the driving reservoirs j = 1,2, ...,j) 
to the reference reservoirs fc = 1,2, ...,d. The probabilities p(n) = {pa-{n)} that the subsystem has evolved to some 
coarse-grained state a while n electrons have been transferred during the time interval t are ruled by the master 
equation: 

9tPt(n) = L-pt(n) (124) 

with the operator 



L-EfL 



(+) --2- 



LinLi-^-- 



(125) 



where L„ is the matrix with the rates of the transitions Up —^ Up zt 1, while L„ is the matrix with the rates of the 
other transitions and all th e loss rates. Examples of such transport processes in quantum dots have been studied for 
instance in Refs. |l04l . ll05| . The cumulant generating function is defined as 

Q(r7)= lim -- In (e"^"), (126) 

where the statistical average is carried out over the probability distribution pt(n). As a consequence, the generating 
function turns out to be given by the eigenvalue problem 

L(»j)-v=-Q(r?)v (127) 

for the modified matricial operator 

L(7,) = e'"-" L e+''-" - Y. ["-i^^ e-"- + if^ + LJ,"^ e+"-] . (128) 

p 

Such matricial operators may obey the following symmetry relation 

M"^-L(7?)- M = L(A-?7)'^ (129) 

where "^ denotes the matricial transpose and M is the matrix of the thermal distribution for the subsystem at 



equilibrium with the d reference reservoirs [28|, llOSj . If this symmetry relation holds for the operator (|128p . this 
property extends to all its eigenvalues and, in particular, to the leading eigenvalue which gives the cumulant generating 
function. In this way, the current fluctuation theorem is demonstrated for the subsystem in the nonequilibrium steady 
state specified by the affinities A: 

Qir,) = QiA-rj). (130) 

Accordingly, the results of Section |V] apply here also for the response coefficients. 

An equivalent expression of the current fluctuation theorem is given in terms of the probability Pt(n) that n 
electrons have been transferred during the time interval t between the reservoirs under the stationary conditions A. 
Equation (|130p implies that opposite fluctuations of these numbers have probabilities obeying the symmetry relation: 



Pt{-n) - 



for t^oo . (131) 



At equilibrium where the affinities are vanishing A = 0, the probabilities of opposite fluctuations are equal so that we 
recover the principle of detailed balancing. In contrast, a directionality manifests itself out of equilibrium where the 
non- vanishing affinities introduce a bias between the probabilities. Soon, one of the probabilities becomes dominant 
as time increases and the currents tend to flow in one direction so that the mean currents 

J = ^(0)=limMl (132) 

are no longer vanishing. 



20 

B. Thermodynamic entropy production 

A consequence of the current fluctuation theorem is the non- negativit y of the entropy production, which is known to 
be eq ual to the sum of the affinities multiplied by the mean currents |69l [74l - [76| . By Jensen's inequality (e^"^) > e^*-^^ 
[T06J . we find that 

(e"^-"), > c-<^-">' (133) 

Consequently, the entropy production is non negative because 

1 ^ = A • J = hm i A • {n)t > lim -1 In (e-^") , = Q(A) = Q(0) = . (134) 

The last equality is the consequence of the definition (|126p for the cumulant generating function at »/ = 0. The second 
law of thermodynamics can thus be deduced from the current fluctuation theorem. 

The entropy production vanishes at equilibrium and is positive out of equilibrium where energy is dissipated. 
In order to drive a particular current with the other currents, energy should thus be supplied and the second law 
constitutes a limit to the efficiency of energy transduction. Powering the process 7 by the other ones requires that the 
corresponding mean current is opposite to its associated affinity: A^J^ < 0. In this case, a thermodynamic efficiency 
can be defined and the second law implies that it cannot reach values larger than unity: 

0<V.^- ^'^''^' , <1- (135) 



Such th ermo dynamic efficiencies have been considered for molecular motors |107l | as well as for mass separation by 
effusion f ioi l. 



C. The case of effusion processes 

The current fluctuation theorem has also been established for effusion processes jlOSl Il09l | . Effusion processes are 
the classical analogues of electron quantum transport through a constriction, such as the quantum point contact of 
Fig. [3K- In effusion, two gases at different pressures and temperatures are separated by a wall with a pore which 
is smaller than the mean free path so that the flow of particles across the pore is essent ially ballistic. The master 
equation of effusion has been known in classical kinetic theory for a century [llOl llllj . The particles of energy 
£ — mw^ j1 are crossing a pore of cross-section area a with the rate: 



a n.i 



W,{e)^^^l3^'ee-^'' for j=L,R (136) 

VSvrm 

if they come from the left- or the right-hand reservoir at the inverse temperature j3j = (fceTj)"^ and the particle 
density rij related to the pressure by Pj = rijk^Tj. These rates satisfy the condition 



W^ie) 



geAE+A« (^37^ 



in terms of the affinities (|104p and (|105p . In effusion, the cumulant generating function of the energy and particle 
currents is given by 

Q(e, r,) = / de [Wi^ie) (l - e-^?-") + W^{e) (l - e^«+'')] . (138) 



in terms of the transition rates ()136p jlOSl Il09| . By their property (I137|) . the generating function obeys the same 
symmetry relation (jl03p as in the quantum case. 

The generating functions (I102p and (|138p of the quantum and classical processes can be compared. At a given 
energy e for ^ = 0, the generating function has the form (|108l) for the quantum process. In the limit where p± <^ 1, 
it becomes 

Q{v) - rp+ (1 - e-") + rp_ (1 - e") (139) 

which is the classical form appearing in Eq. (I138P with the identification rp± — VFl,r(£). The expression (I139P is 
characteristic of the combination of two independent Poisson processes, i.e., the particles can be transferred in both 
directions. 
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D. Statistics of histories and time reversal 

Another time-reversal symmetry relationship concerns the statistical properties of the histories or paths followed 
by a system under stroboscopic observations at some sampling time At. Such observations generate sequences of 
coarse-grained states such as 

OJ = UJ1LU2 ■ ■ -UJn (140) 

corresponding to the successive times tj = j At with j — 1,2, ..., n. This history or path has a certain probability P{(x>) 
to happen if the system is in the stationary state corresponding to the affinities A. Because of the randomness of the 
molecular fluctuations, these path probabilities typically decrease exponentially at some rate h that characterizes the 
temporal disorder in the process. This characterization applies to stochastic processes as well as chaotic dynamical 
systems, for which the temporal disorder h is called the Kolmogorov-Sinai or dynamical entropy (tI. H, llSl , lll2J . 
In nonequilibrium stationary states, the time-reversed path 

W^ = a;„ • • • a;2Wi (141) 

is expected to happen with a different probability decreasing at the different rate h^ now characterizing the temporal 
disorder of the time- reversed paths [33|- The remarkable result is that the difference between the disorders of the 
time-reversed and typical paths is equal to the thermodynamic entropy production '37'| 

^M = /,R_/i>0 (142) 

Ub dt - ^ ' 

The second law of thermodyn amic s is satisfied because this difference is known in mathematics to be a relative entropy, 
which is always non negative |106{ . At equilibrium, detailed balancing holds so that every history and its time reversal 
are equiprobable, their temporal disorders are equal h = h^ , and the entropy production vanishes. This is no longer 
the case away from equilibrium where the typical paths are more probable than their time reversals. Consequently, 
the time-reversal symmetry is broken at the level of the statistical description in terms of the probability distribution 
of the nonequilibrium stationary state. In this regard, the entropy production is a measure of the time asymmetry in 
the temporal disorders of the typical histories and their time reversals. As a corollary of the second law, the typical 
histories are more ordered in time than their corresponding time reversals in the sense that h < h^ in nonequilibrium 
stationary states |40l |. 

In the case of effusion processes, the temporal disorders of the histories and their time reversals are characterized 

by 



h'' = 



e \ f°° r°° 

In^^)/ de[WL{e) + Wu.{e)]- I de[WL{e)\nWL{e) + WR{e)\nWR{e)] + 0{At) , (143) 
O'^AiAt)/ de[WL{e) + Wu.{e)]- I de[WL{e)\nWR{e) + WR{e)\nWL{e)] + 0{At) , (144) 



in terms of the transition rates (jl36p . the sampling time At, and the coarse graining in energy Ae. The first term 
of these quantities is a feature of stochastic processes that are continuous in time and energy. For such processes, 
randomness manifests itself on arbitrarily small time scales At and energy scales Ae so that the temporal disorders 
increase as AeAi — >■ 0. 

The difference between these two quantities is equal to the entropy production because 



h'' -h = 



h= r de [lyL(e) - W^R(e)] In ^44 = Ae Je + A^ Jn = ^'^ . (145) 

Jq WR,(e) fcB dt 

in the limit where AeAt -^ 0. 

The validity of the formula (|142|) has also been verified in experiments where the nonequilibrium constraints are 
imposed by fixing the currents instead of the affinities [4l|, |42] . 

Similar considerations have been developed for quantum systems |43l . Ill3l | . We n otice that quantum mechanics 
naturally limits the randomness on the scale where AeAt — 2TrH, as explained in Ref. J113J |. 

VIII. TRANSPORT IN ELECTRONIC CIRCUITS 

Several types of electronic circuits are considered such as single-electron transistors, quantum dots with resonant 
levels, or quantum dots capacitively coupled with a quantum point contact |3l. \4. l90l| . 



22 

A. Quantum dot with one resonant level 

A quantum dot is a kind of artificial atom with quantized electronic levels. Because of the tunneling to the reservoirs, 
these levels are broadened into resonances characterized by a lifetime. We suppose that the transport process involves 
a single resonant level, which may thus be either empty or occupied. In this case, the process is ruled by the master 
equation (|124p with a modified operator (II28P given by 

■-W^ r...-r,^^^ u _u (146) 
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in terms of the following charging and discharging rates: 

ap^Tpfp and 6p == Tp (1 - /p) for p = L,R. (147) 

The rate constants Fp are proportional to the square of the interaction parameters between the quantum dot and the 
reservoir p, as well as to the density of states of the reservoir p. The Fermi-Dirac distributions are given by 

/o = TTr — rrrr r for p = L,R (148) 

where £q is the bare energy of t he do t level, AC/q is the electrostatic charging energy, and fip = eVp arc the chemical 
potentials of the reservoirs 0, Ij, lll4J [ . 

The cumulant generating function is given by the eigenvalue problem (I127P for the matrix (jl46p as 
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(149) 



This function obeys the symmetry relation Q{r]) = Q{A — rj) whereupon the current fluctuation theorem (|13ip is 
satisfied with the affinity [7^ 

^^ln^=/3(^L-MR)=/3eV^- (150) 

aROL 

B. Capacitively coupled circuits 

In order to perform the full counting statistics of electron transport in quantum dots, the main circuit should 
be monitored by a secondary circuit which is made of a quantum point contact (QPC) [59146 1|. Both circuits are 
capacitively coupled so that the current in the QPC is sensitive to the occupancy of the QDs. This latter modulates 
the QPC current by Coulomb repulsion. In this way, the quantum jumps in the occupancy of the QDs can be 
experimentally observed. Typically, the electric current in the QPC is seven or eight orders of magnitude larger than 
the current in the circuit with the QDs (see Table U). From the viewpoint of quantum measurement, the QPC plays 
the role of a measuring device and the QDs the small observed quantum system. Because of the large ratio between 
their currents, the measuring device is essentially in a classical regime. Although the ratio of the noise to the current 
is negligible in this regime, the noise gets larger as the current increases. As a consequence, the QDs themselves are 
subjected to an important noise, which affects their kinetics. In particular, the charging and discharging rates of the 
QDs are modified by the presence of the large current in the QPC. Therefore, the Coulomb interaction between both 
circuits also causes a back action of the QPC onto the QDs. 

Several types of circuits have been considered. In Ref. [S^, a single QD with a strong bias voltage is monitored 
by a QPC. The QD is successively empty or occupied so that the QPC current is observed to jump between both 
corresponding values. This setup allows to carry out the full counting statistics of the occupancy of the QD. Since 
the bias voltage is large fin, ^ Mir, the transition rates of the backward fluctuations can be assumed to be negligible. 
Under these conditions, the statistics of the fiuctuating current can be inferred from the statistics of the occupancy. 
However, the backward fiuctuations of the current are also negligible in this fully irreversible regime. Therefore, the 
current fluctuation theorem cannot be tested for lack of backward fluctuations. 

In order to overcome such limitations, a setup with two QDs has been considered in Ref. [6l|. The two QDs are 
positioned at different distances from the QPC so that different Coulomb interactions U and U' exist between the 
QPC and both QDs, as shown in Fig. [SJd- If an electron travels from the left- to the right-hand QD, the Coulomb 
repulsion is successively larger and smaller, or vice versa for the reversed motion of an electron. Therefore, this setup 
is suitable for bidirectional counting statistics. Remarkably, the symmetry relation Pt{n) = Pi(— n) exp(AQD ti) is 
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FIG. 5: (a) Schematic representation of a single quantum dot (QD) monitored by a quantum point contact (QPC) to which 
it is capacitively coupled by the Coulomb interaction U, as in the circuit of Ref. [5y|. (b) Schematic representation of two 
QDs monitored by a QPC asymmetrically positioned with respect to the QDs in order to get different capacitive couplings of 
the QPC to the two QDs, as in the circuits of Refs. [6ll - [6^ . The Coulomb repulsion U is larger than [/', which allows the 
bidirectional counting of electrons in the circuit with the two QDs. 





single QD [59] 


double QD [61J 


T 


350 mK 


130 mK 


Iqd 


1.3 X 10"^" A 


6.5 X lO-i'^ A 


Vqu 


2.7 mV 
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Hqd 


3.4 X 10"^^ W 


2.0 X 10-2° W 


>^QD 


792 Hz 


406 Hz 


^QD 


90 


27 


-^QPC 


4.5 X 10"^ A 


1.2 X 10-** A 


Vqpc 


0.5 mV 


0.8 mV 


nqpc 


2.3 X 10-12 W 


9.6 X 10-12 ^ 


Jqpc 


2.8 X 10"' Hz 


7.5 X lO^" Hz 


Aqpc 


17 


71 


Iqpc 
Iqd 


3.5 X 10^ 


1.8 X 10** 



TABLE I: Electronic temperature T, mean electric currents la, voltages Va, dissipated powers Ila = Vala, mean electron 
currents J^ = 7a/|e|, affinities A^ ~ \e\Va/{kBT) in the quantum dot (a = QD) and the quantum point contact {a = QPC), as 
well as the ratio of QPC to QD currents, for the experiments reported in Ref. [53] with a single QD and Ref. 61] with a double 
QD. According to Eq. (|134p . the thermodynamic entropy production in the whole system is equal to diS/dt = fca "^^ Ac J a = 
X^Q n^/r, which is dominated by the energy dissipation in the QPC under these experimental conditions. 



observed but, depending on the experimental conditions JGl'-'GSj, the effective affinity Aqo may significantly differ from 
the thermodynamic affinity ^qd = \e\VQu/{kBT) driving the current in the QD. In Refs. ^,£2,], the gate voltages 
of the QDs are tuned so that transitions occur between charge states where the QDs are singly or doubly occupied 
and the effective affinity Aqd — 3 differs by one order of magnitude from the thermodynamic affinity Aqd = 27 (see 
Table |T]) . In Ref. |63j , the experiment has been carried out closer to equihbrium with an optimized sample design 
and gate voltages tuned so that transitions occur between charge states where the QDs are empty or singly occupied. 
For the bias voltage Vqd = 20 ^V and the electronic temperature T = 330 mK, the thermodynamic affinity is equal 
to Aqd = 0.7 while the effective affinity is observed to be Aqu — 0.5. The differences observed in the experiments 
reported in Refs. |6ll - [63| arise because the QPC may induce back-action effects onto the dynamics of the QDs under 
some conditions, while the limited bandwidth of the charge detection also affects the counting statistics. 

In order to understand the back-action effects, the capacitive coupling between the QDs and the QPC should be 
taken into account. For th is purpose, th e simple model of Fig. [5^ with two capacitively coupled QDs in two parallel 
circuits may be considered |l04ll05l . lll6{ . The system is composed of d = 2 disconnected circuits and r = 4 reservoirs, 
so that the nonequilibrium steady states are specified with p = r — d ~ 2 affinities: 

Ai = ^(^iL - Mir) and A2 = fi{^J.2l. - ^J■2R) ■ (151) 

The system can be described in terms of the probabilities Piyj^i,^ (f^l , n2) that each quantum dot (a = 1 or 2) is occupied 
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FIG. 6: (a) Schematic representation of two quantum dots in parallel. Each quantum dot is coupled to two reservoirs of 
electrons. Moreover, both quantum dots influence each other by the Coulomb electrostatic interaction U . (b) Simulation with 
Gillespie's algorithm yj.5] of the detector current in circuit No. 2 measuring the occupancy of the QD No. 1. The parameter 
values are given by /3ei = 0, ^£2 = 35, PU = 32.8, /3/iiL = 25, ^^iR = 0, Til = Tm = f il = f m = 1, ^^2l = 70, P^l2n = 0, 
and r2L = r2R, — r2L = r2R — lO*. The affinities of both circuits are Ai = 25 and A2 = 70. The effective affinity of the circuit 
No. f is y4i = f .f7. The mean value of the QD current is J\ ~ 0.17 electrons per unit time. The mean value of the secondary 
current is J2 ~ 4.8 x 10^ electrons per unit time. The QD is empty (respectively occupied) when the secondary current takes 
the value 5 x 10^ (respectively 4.5 x 10^). Adapted from Ref. ji)5|]. 



by i'q = or 1 electron, while Ua electrons have been transferred from the left- to the right-hand reservoirs in the 
circuit a during the time interval t. The stochastic process is ruled by the master equation (|124p with the operator 
p25p . For a = 1, 2 and p = L, R, the charging and discharging rates are given by 



either 



or 



T^apfap and bap = Tapil - lap) with faf. 



1 



O'o.p — t apja 



pjap 



and 



bap = Tap{l - fap) with fap = 



g/3(£Q-Map) + 1 ' 

1 



(152) 
(153) 



if the other quantum dot is eith er empty or occupied. The effect of the tunneling barri er ca pacitances is neglected in 
this simplified model jl04 Il05| . We notice that their effect has been analyzed in Ref. [116J | . 

Now, the cumulant generating function (I126P is obtained by solving the eigenvalue problem (I127P and it obeys the 
current fluctuation theorem (I130p for the two currents flowing across the parallel circuits: 



Qim,m) = 2(^1 -'717^2 - m) 



(154) 



with respect to the two affinities ()15ip . However, this result does not imply that the symmetry relation holds for the 
single current flowing in one of the circuits. The generating function of this single current is given by Q(j]i,0) and 
the symmetry does not hold in general, Q(?7i, 0) 7^ Q{Ai — ?yi, 0), as shown in the counterexample of Fig. [7^. 
Nevertheless, under some circumstances, the single-current fluctuation theorem 



2(771,0) = 2(^1-7,1,0) 



(155) 



holds with respect to an effective affinity Ai, which can be calculated. 

This symmetry holds with the effective affinity equal to the thermodynamic affinity, Ai = Ai, in a few cases: (1) if 
the other circuit is at equilibrium, i.e., if A2 = 0; (2) if there is no Coulomb interaction between both circuits, t/ = 0, 
in which case the two circuits are independent of each other; (3) in the limit where the Coulomb interaction is very 
large, U = c». Indeed, the probability that both quantum dots are occup ied is vanishing in this limit, which reduces 
the Markov process to three states and implies the symmetry (|155[) (l05J | . 

The single-current fluctuation theorem (|155|) also holds in the limit where the rate constants of one circuit are much 
larger than in the other circuit: 



^2p,^2p ^ Fip 



,r 



ip 



for 



p = L,R 



(156) 



as in the experiments of Refs. f59|,|61| where the ratio between the QPC and QD currents is many orders of magnitude 
larger than unity (see Tabic H]). In the system of Fig. [5^, the fast QD in the secondary circuit rapidly jumps between 
its two states 7^2 = and 7^2 = 1 so that the charging and discharging rates of the slow QD are effectively averaged 
over both states. The conditional probabilities Pj^^l^i ^^at the second QD has the occupancy 7^2 provided that the 
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first QD is in the state vi are given by 



p _ &2L + b2R 

a2L + a2R + ^2L + &2R 



^o|i — 
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Pill- 



a2L + a2R 



a2L + a2R + &2L + &2R 

a2L + a2R 



a2L + a2R + 62L + ^2R a2L + a2R + &2L + &2R 

Therefore, the effective charging and discharging rates of the first QD are given by the following averages: 



flp — aipFo|o 



aipPi|o 



and 



bp — foip^oii + ^ip All 



with 



p = L,R. 



(157) 
(158) 

(159) 



Under the conditions (I156p . the Markovian process is similar to the one ruled by the operator (|146p but with these 
effective charging and discharging rates. Accordingly, the cumulant generating function has the form (jl49p with the 
rates (|159p so that the single-current fluctuation theorem (|155p is indeed satisfied with the effective affinity: 



Ai = ln 



aL^R 
aR^L 



(160) 
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FIG. 7: Tlie cumulant generating function Q{rii , 772 = 0) for the circuit of Fig. |6^ versus the counting parameter r/i : (a) together 
with the symmetric function with respect to the effective affinity Ai = 16.8356 (dotted-dashed line) for /3[/ — 10 and T2h ~ 
T2R = r2L = r2R = 2; (b) for different values of the electrostatic coupling parameter f5U and r2 = r2L ~ r2R — T2h ~ r2R ~ 
100. In both cases, the other parameters take the following values: /?ei = 10, /3e2 ~ 35, /3/iiL = 25, /3^ir = 0, /3^2l = 70, 
/3/i2R ~ 0, and Fi = Fil = Fir = Fil = Fir = 1. In these plots, the effective affinity is given by the second root of the func tion 
at rji = Ai ^ 0. The thermodynamic affinity of the secondary circuit is A2 — /?(/i2L — M2r) ~ 70. Adapted from Ref. [l05l |. 

Figurellh) shows not only that the single-current generating function Q(?7i, 0) is nearly symmetric already if r2/ri = 
100, but also that the symmetry holds with respect to an effective affinity Ai which varies with the electrostatic 
interaction l3U . If the Coulomb interaction vanishes l3U = 0, both circuits are decoupled so that the thermodynamic 
affinity Ai — 25 is recovered. However, the effective affinity decreases if the electrostatic interaction increases till 
/3U = 20, as seen in Fig.[7j3. For larger interaction /3U, the effective affinity (|160p recovers its thermodynamic value 
in the limit U —>■ 00, as explained here above ^105i] . 

The conclusion is that the secondary circuit monitoring the quantum jump dynamics in quantum dots may have 
important back-action effects on the full counting statistics. In general, the current fluctuation theorem holds at 
the fundamental level of description for all the currents interacting in the system. In this regard, the secondary 
circuit performing the measurement of the quantum state in the first circuit is part of the total system and cannot 
be separated in general. Nevertheless, the symmetry relation of the fluctuation theorem is still observed to hold for 
a single fluctuating current, but with an effective affinity that may differ from the thermodynamic value due to the 
back-action effect of the monitoring circuit onto the observed circuit. The effective affinity (|160p may drop down to 
a value more than one order of magnitude lower than the thermodynamic affinity [ 105] . Similar conclusions have 
been reached with related approac hes where the averages (|159p of the transition rates have been carried out in the 
framework of the P(£;)-theory [6^. [TTtI . [TTsI . 

A general remark is that quantum measurement implies the dissipation of energy. In order to observe the quantum 
jumps of the QD with a short enough time resolution, the current and the dissipated power in the monitoring circuit 
should be higher than in the QD. If the QD state is monitored with a sampling time Ai, the secondary circuit playing 



the role of the detector should have a current J2 > (At)~ . Hence 

II2 = kBTA2J2 ^ kBTA2{At)~^. The higher the time resolution, the higher the dissipation rate. 

experimental conditions are variously influenced by the bandwidth of charge detection. 



the dissipated power should be bounded by 

Therefore, the 
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C. Coherent quantum conductor 

The current fluctuation theorem also applies to transport in coherent quantum conductors such as Aharonov-Bohm 
rings in a magnetic field @t3- We suppose that the circuit is connected to two reservoirs so that there is a single 
affinity. In this case, the results of Section |V] are the following. The particle current is expanded in powers of the 
affinity according to Eq. (|65l) : 

JiA;B) = L{B)A+^M{B)A'^ + --- (161) 

Similar expansions can be introduced for the diffusivity (j63p and the third cumulant (j64p : 

D{A; B) = DoiB) + Di{B) A + --- , (162) 

C{A;B)^CoiB) + --- . (163) 

As shown in Scction|Vl all these coefficients are interconnected by the following relations. For linear response and the 
second cumulant, Eqs. ([68|) - (|69)) give 



LiB) = L{^B) = Do{B) . (164) 

For the third cumulant and the second response coefficient, Eqs. ((7^ - (fM)) give 

Co{B) = -Co{-B) = 2 [Di{B) - Di{-B)] , (165) 

M{B) + M{-B) ^ 2 [Di (B) + Di{-B)] , (166) 

M{B)^2DiiB)-^Co{B). (167) 

Similar relations hold at higher orders |52| . 

In order to compare with usual quantities, we notice that the electric current is equal to I{V; B) = eJ{A; B) and 
the affinity A is related to the voltage V by 

eV 
kbT 

Moreover, the noise power ()114|) at zero frequency S — S{uj = 0) is proportional to the diffusivity according to 
Eq. (|115p . It is standard to expand the electric current and the noise power as follows: 

I{V;B) = G,{B)V +^G2{B)V^ + ■ ■ ■ , (169) 

SiV-B) = SoiB) + SiiB)V + --- . (170) 

Comparing with the expansions (I16ip and (|162l) . we get 

e'^L^keTGi, e' M ^ {kBTf G2 , ... (171) 

Ae^Df)^So, Ae^ Di ^ keT Si , ... (172) 

Accordingly, we recover the Johnson-Nyquist and Casimir-Onsager relations [93| 

So{B) = 'ikBTGi{B) and Gi{B) = Gi{-B) . (173) 

For the nonlinear response properties, the following relations can be deduced 

SiiB) + Sii-B) = 2kBT[G2{B) + G2{-B)] , (174) 

SiiB)-Sii-B) = 6kBT[G2{B)-G2{-B)] , (175) 

at the level of the third cumulant Co{B), as well as higher-order relations. 

The expression (|174[) is only the consequence of the second equality in Eq. ((58)) . Q{A, A; B) — 0, which has been called 
the 'global detailed balancing condition' and which is weaker than the current ffuctuation theorem itself [1191] . The 
expression (|175p is the consequence of the current fluctuation theorem and, thus, of the underlying microreversibility 
in the sense that QH{B)Q~^ = H{—B). We notice that the quantity (|175p is proportional to the magnetic asymmetry 
given by the third cumulant (jl65p . 
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The symmetry relations (|174[) and (1175 1 ) h ave been tested in recent experiments on a circuit with an Aharonov- 
Bohm ring in an external magnetic field [6J, [6a| . The circuit is fabricated by local oxidation on a GaAs/ AlGaAs 
2DEG and the ring has a diameter of about 500 nm. The electronic temperature T = 125 mK is determined with 
the Johnson-Nyquist fluctuation-dissipation relation. The current is observed as a function of the applied voltage 
V, the gate voltage Vg, and the external magnetic field B. The noise power is measured with the cross-correlation 
technique between the signals from two sets of resonant circuit and amplifier, which are external to the circuit. Beyond 
linear response, the coefficients Si and G2 have been measured and they have complex dependences on the magnetic 
field and the gate voltage. Remarkably, the proportionalities predicted by Eqs. (|174p and (I175P are confirmed. The 
proportionality constants are observed to be 12.00 ± 1.96 instead of 2 in Eq. (|174|) and 9.66 ± 1.32 instead of 6 in 
Eq. (|175p [65| . The antisymmetric relation (I175P is in better agreement with theory than the symmetric relation (11741) . 
The reason for the considerable deviation in the case of the symmetric relation (I174p is not explained in Refs. [64, 65]. 
The amplifiers used to measure the noise power may have back-action effects. 

The fundamental questions associated with these issues deserve further experiments to understand the origins of 
the quantitative deviations and to test predictions at the higher orders if the concerned quantities are experimentally 
accessible. 

IX. CONCLUSIONS 

This chapter has been devoted to recent advances in the nonequilibrium statistical mechanics of quantum systems. 
For more than a century, statistical mechanics has offered the conceptual framework to formulate the various statistical 
aspects of mechanics. These aspects play a prominent role in quantum mechanics, which is an initial-condition 
theory as well as classical mechanics. In such theories, the non-reproducibility of random phenomena is naturally 
explained by the arbitrariness of initial conditions in due respect to the principle of causality. This is the case in open 
systems under nonequilibrium conditions sustaining currents of energy and particles between reservoirs at different 
temperatures and chemical potentials. After transients, the system evolves towards a nonequilibrium steady state 
which describes the different properties such as the mean currents and their statistical cumulants. The set of all these 
cumulants characterize the full counting statistics and, thus, the large-deviation properties of the fluctuating currents, 
which have been the focus of recent advances in nonequilibrium statistical mechanics. 

These advances find their origins in earlier work on chaotic dynamical systems modeling transport properties such 
as diffusion, viscosity or heat conductivity |8l4lq. In this context, large-deviation properties have been studied for 
physical quantities fiuctuating in time. Previously, the large-deviation properties were mainly considered in equilibrium 
statistical mechanics to infer thermodynamics from the spatial fluctuations of observables in the large-size limit. In 
dynamical systems theory, large-deviation properties have been introduced in close relation with the concepts of 
Lyapunov exponent, Kolmogorov-Sinai entropy per unit time, and fractal dimensions [7|. Furthermore, relationships 
have been established between these quantities characterizing dynamical chaos and the transport properties |8l4l8l| . In 
this context, the consideration of time-reversal symmetry for nonequilibrium steady states has opened new perspectives 
for a fundamental understanding of nonequilibrium systems at small scales where the microscopic degrees of freedom 
manifest themselves as fluctuations. Several types of time-reversal symmetry relations have been established |2l| - 
139( 1 . On the one hand, so-called fluctuation theorems have been proved which compare the probabilities of opposite 
fluctuations for nonequilibrium quantities of interest. On the other hand, the thermodynamic entropy production 
has been shown to result from the breaking of time-reversal symmetry at the statistical level of description. These 
advances have turned out to play a fundamental role in the understanding of many nonequilibrium nanosystems such 
as molecular machines. 

More recently, time-reversal symmetry relations have been extended to quantum systems and, in particular, to 
electron quantum transport in semiconducting nanodevices |43l - [58| . Several approaches have been followed. 

Functional relations have been established which provide a unifying framework to deduce many results such as the 
Casimir-Onsager reciprocity relations, the Kubo formula of linear response theory, the quantum version of Jarzynski's 
nonequilibrium work equality, as well as their generalization beyond linear response or to novel physical contexts [5l| . 

Moreover, a theoretical approach has been developed for the transport of energy and particles across a driven 
open quantum system between reservoirs under nonequilibrium conditions. In quantum mechanics, the issue of 
measurement is always subtle because of the non-separability of a system into its different parts. Here, quantum 
measurements at the initial and final times are required in order to determine the amounts of energy and particles 
transferred between the reservoirs during the lapse of time when the system is driven by time-dependent forces. Using 
the time-reversal symmetry of the microscopic dynamics or its extension in the presence of an external magnetic field, 
symmetry relations have been obtained which are quantum versions of the transitory current fluctuation theorem. 

In the long-time limit, this theorem leads to the stationary current fluctuation theorem. An essential aspect is that, 
for nonequilibrium stationary conditions, the symmetry relation should hold with respect to the affinities or thermo- 
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dynamic forces due to the differences of temperatures or chemical potentials between the reservoirs. Remarkably, a 
proposition guarantees that the cumulant generating function of the fluctuating currents indeed depends only on the 
differences between the counting parameters associated with the reservoirs between which the thermodynamic forces 
are maintained J52| . Thanks to this proposition enunciated in Section IIVI the stationary current fluctuation theorem 
can be established for all the currents flowing across open quantum systems. 

Close to equilibrium, known results of linear response theory such as the Casimir-Onsager reciprocity relations and 
fluctuation-dissipation formulae can be deduced from the stationary current fluctuation theorem. The fact is that 
this theorem also holds arbitrarily far away from equilibrium, allowing the generalization of these results to nonlinear 
response properties [5^ [77H79J . 

These advances concern in particular electron quantum transport in mesoscopic semiconducting circuits. At low 
temperature, electrons have ballistic motions on the size of such circuits, which form electronic waveguides. Many 
transport properties of mesoscopic conductors can thus be understood by supposing that the electrons are independent 
of each other and undergo scattering in the effective potential formed by the semiconducting circuit. In this approx- 
imation where electron-electron and electron-phonon interactions are neglected, the cumulant generating function of 
the aforementioned theorem can be expressed with the Levitov-Lesovik formula in terms of the scattering matrix of 
the independent electrons on the obstacles of the circuit. Therefore, this generating function has the time-reversal 
symmetry of the fluctuation theorem with respect to the affinities driving the currents across the circuit. This fun- 
damental result has many consequences, in particular, concerning the full counting statistics of electrons and the 
relationships between conductance and noise power beyond linear response. Recently, these consequences have been 
investigated experimentally |62| - [65| . Since the full counting statistics is performed by the capacitive coupling to a 
secondary circuit, back-action effects should be taken into account to interpret the experimental results. Accordingly, 
the symmetry relation for the single observed current holds with respect to an effe ctive affinity which differs from the 
thermodynamic affinity because of the noise coming from the secondary circuit |62l . ll05| . In this context, new opportu- 
niti es ar ise for a better understanding of the intriguing fact that quantum jumps can be observed continuously in time 
99l - llQl| . For electronic conduction in an external magnetic field, the current fluctuation theorem predicts relations 



beyond linear response between conductance and noise power [50l |52| | . Experimental observation on Aharonov-Bohm 
rings designed with GaAs semiconductor are in semiquantitative agreement with the theoretical predictions, possibly 
because of back action from the noise measurement [6J, |65| . 

In conclusion, the advent of time-reversal symmetry relations among the dynamical large-deviation properties has 
opened new perspectives in our fundamental understanding of nonequilibrium quantum systems. These relations can 
be extended to many nonequilibrium systems beyond electron quantum transport. In particular, they may apply 
to fermion and boson quantum transport in the physics of cold atoms or molecules or in quantum optics. The 
properties which can be deduced from the symmetry relations notably concern the properties of energy transduction 
when several currents are coupled together in thermoelectric, photoelectric, or mass separation devices. In a broade r 
perspective, time-reversal symmetry relations are also considered in relativistic and gravitational systems |l20lll21j . 
Similar relatio ns ca n be envisaged for other discrete symmetries or their combinations, e.g., in the context of quantum 
field theory fi\.\l2^. 
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